Is the current lack of solar activity only skin deep? 
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ABSTRACT 

The Sun is a variable star whose magnetic activity and total irradiance vary on a timescale of ap- 
proximately 1 1 years. The current activity minimum has attracted considerable interest because of its 
unusual duration and depth. This raises the question: what might be happening beneath the surface where 
the magnetic activity ultimately originates? The surface activity can be linked to the conditions in the 
solar interior by the observation and analysis of the frequencies of the Sun's natural seismic modes of 
oscillation - the p modes. These seismic frequencies respond to changes in activity and are probes of con- 
ditions within the Sun. The Birmingham Solar-Oscillations Network (BiSON) has made measurements 
of p-mode frequencies over the last three solar activity cycles, and so is in a unique position to explore the 
current unusual and extended solar minimum. We show that the BiSON data reveal significant variations 
of the p-mode frequencies during the current minimum. This is in marked contrast to the surface activity 
observations, which show little variation over the same period. The level of the minimum is significantly 
deeper in the p-mode frequencies than in the surface observations. We observe a quasi-biennial signal in 
the p-mode frequencies, which has not previously been observed at mid- and low-activity levels. The stark 
differences in the behavior of the frequencies and the surface activity measures point to activity-related 
processes occurring in the solar interior, which are yet to reach the surface, where they may be attenuated. 

Subject headings: methods: data analysis - Sun: activity - Sun: helioseismology - Sun: oscillations 



1. Introduction 

The level of the Sun's magnetic activity is observed 
to vary on an 1 1 -year time scale and we are currently in 
the minimum between cycles 23 and 24. The current 
solar minimum is attracting a great deal of attention 
as it is proving to be quite unusual. Observations of 
surface and atmospheric effects, such as the number of 
visible sunspots, the rate of occurrence of solar flares 
and the strength of the solar wind, highlight just how 
quiet the Sun is^]. 

The Sun's activity cycle influences everyday life on 
the Earth. The rate of occurrence of solar flares is de- 
pendent on the number of spots on the Sun's surface, 
and large solar flares can disrupt satellite communica- 
tions and cause power outages. Coronal mass ejections 



(CMEs), which are another source of radiation that can 
disrupt life on the Earth, are still being observed regu- 
larly on the Sun, even in this unusual solar minimum 
(see the STEREO COR1 CME catalog^). Cosmic rays, 
which are a significant space radiation hazard, are an- 
ticorrelated with the solar cycle. The level of solar ac- 
tivity is possibl y correlated to the Earth's c limate (see, 
for example, lLockwood & Frohlichl 120071 and refer- 
ences therein). Space-weather groups use predictions 
of solar cycles to anticipate the amount of orbital drag 
experienced by satellites. 

The next solar cycle has already proven difficult to 
predict as the current solar mi nimum i s lasti ng signif- 
icantly longer than expected. IPesnelll ( 2008 ) reviews 
50 predictions for the size and timing of cycle 24 and 
finds a wide range of results, especially in compari- 



j http://science.nasa.gov/headlines/y2009/01apr.deep S olarminimum.htm ^ |http://corl.gsfc.nasa.gov/docs/prelim_ev e nts/CORl P reUmCMErate _Feb2Q09.pdf 
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son to prediction s made before the previous solar cy- 
cle (IJoselyn et al.lll997l) . For example, predictions of 
when cycle 24 will reach its maximu m range from 
2009 D ecember dMaris & Oncicall2006l) to 2014 De- 
cember dTsirulnik. Kuznetsova & Oraevskv 1997 ). In 
fact, the official NOAA, NASA and ISES Solar Cycle 
24 Predict ion Panel, wh ich studied the predictions col- 
lated by |Pesnei3 ( 2008 ). failed to reach a consensus on 
when the peak of cycle 24 will occur and how active 
the upcoming cycle will be. Meanwhile, the number 
of predictions for cycle 24 is ever increasing. 

Surface measures of the Sun's activity, such as the 
number of sunspots, which are used to aid cycle pre- 
dictions, indicate that we are still in an extended solar- 
cycle minimum. We ask the question: can we learn 
anything about this unusual solar minimum from the 
Sun's interior? 

To answer this question we investigate the variation 
with the solar cycle of the frequencies of the Sun's 
natural resonant oscillations, which are known as p 
modes. Solar p modes are trapped in cavities below the 
surface of the Sun and their frequencies are sensitive 
to properties, such as temperature and mean molecular 
weight, of the s olar material. It has been known since 
the mid 1980s dWoodard & Novesl[l985l) that p-mode 
frequencies vary throughout the solar cycle with the 
frequencies being at their largest when the solar activ- 
ity is at its maximum. By examining the changes in 
the observed p-mode frequencies throughout the solar 
cycle, we can learn about solar-cycle-related processes 
that occur beneath the Sun's surface. 

Th e Birmingham Sola r-Oscillations Network (Bi- 
SON: IChaplin et al. 1996) makes Sun-as-a-star (unre- 
solved) Doppler velocity observations, which are sen- 
sitive to the p modes with the largest horizontal scales 
(or the lowest angular degrees, £). Consequently, the 
frequencies measured by BiSON are of the truly global 
modes of the Sun. These modes travel to the Sun's core 
but their dwell time at the surface is longer than at the 
solar core because the sound speed inside the Sun in- 
creases with depth. Consequently, the low-f modes are 
most sensitive to variations in regions of the interior 
that are close to the surface and so are able to give a 
global picture of the influence of near-surface activity. 
BiSON is a network of autonomous ground-based ob- 
servatories that are strategically positioned at various 
longitudes in order to provide as continuous coverage 
as possible of the Sun. BiSON is in a unique position 
to study the changes in oscillation frequencies that ac- 
company the solar cycle as it has now been collecting 



data for over 30 years. However, when the network 
was first established the quality of the data was rela- 
tively poor, in comparison to recent years, because of 
the sporadic coverage of the observations. Here, we 
have been able to analyze the mode frequencies ob- 
served during the last two solar cycles in their entirety. 

2. Analysis and results 

The precision with which p-mode frequencies can 
be determined is directly related to the length of time 
series under consideration. Consequently, p-mode fre- 
quencies are often determined from time series whose 
lengths are of the order of years. However, a com- 
promise must be made here regarding the appropriate 
length of time series for study so that changes in the 
solar cycle can be resolved. Before 1985 April the ob- 
served data are sparse because the early BiSON data 
were collected from just one or two sites and, initially, 
just in summer months in the northern hemisphere. 
Therefore, when considering the data collected be- 
fore 1985 April the p-mode frequencies were obtained 
from seven time series of different lengths, which re- 
flected the availability of data. These data are included 
here for completeness but are not used in the later anal- 
ysis as the data are too sparse to provide reliable re- 
sults. After the third BiSON site, at Carnarvon, West- 
ern Australia, was established in 1985 the duty cycle 
of the BiSON data increased significantly. Therefore, 
after 1985 April 22 the time series were truncated to 
54 days in length, which corresponds to approximately 
twice the rotation period of the solar surface. The anal- 
ysis performed here concentrates on the data obtained 
after this date. 

A standard likelihood maximization method was 
used to fit the powe r spectra of these time series (e.g., 
Chaplin et alJ ll999). enabling the observed mode fre- 
quencies to be determined. We have concentrated on 
the strongest modes of oscillation, which are in the fre- 
quency band from 2100 to 3500 pHz. 

To maintain con sistency with previous work (e.g., 
Chaplin et alj|2007l) a minimum activity reference set 
was determined by averaging the frequencies from 
time series observed during the minimum activity 
epoch at the boundary between cycle 22 and cycle 
23. We then defined the solar-cycle frequency shifts, 
6v„j(t), as the differences between frequencies given 
in the minimum activity reference set and the frequen- 
cies of the corresponding modes observed at different 
epochs and, consequently, different levels of activity. 
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Fig. 1. — Smoothed frequency shifts observed in BiSON data during cycles 21, 22, and 23. The error bars associated 
with the frequency-shift data are due to uncertainties in the fitted oscillation frequencies. The green line is a sinusoid 
which has been fitted to the seismic data. Overplotted in red, and using the right-hand axis, is the F10.7, which has been 
scaled using a linear fit to the observed frequency shifts. 



The size of a frequency shift has a well-known de- 
pendency on frequency and mode inertia and these 
de pendencies were rem oved in the manner described 



Chapl in et al. I J2007h . This allowed a weighted av- 



m 

erage of the frequency shifts observed for each time 
series to be determined, which provided a mean fre- 
quency shift for each epoch, 5v{t). Figure Q] shows 
the frequency-shift data that were observed during cy- 
cles 21, 22, and 23. The data are shown up to 2009 
April. After 1985 April the data plotted in Figure Q] 
have been smoothed over five points, as the mean fre- 
quency shifts were obtained from time series of the 
same length. 

One of the most obvious manifestations of the Sun's 
solar cycle is the variation in the number of sunspots 
on the solar surface, formally recorded as the Interna- 
tional Sunspot Number (ISN). However, the sunspot 
number is not the only variable that can be used as a 
proxy for the solar activity. For example, the radio 
flux emitted from the Sun at a wavelength of 10.7cm 
(hereafter F10.7) also responds to changes in the so- 



lar cycle. Many authors (e.g., lElsworth et al. 1990; 



Howdl2008l and references therein) have determined 



and commented upon the good correlations between 
shifts in the /7-mode frequencies and activity proxies 
such as the F10.7 and the ISN. Therefore, for compari- 
son purposes, the average 10.7 cm flux has been plot- 
ted on top of the observed frequency shifts in Figure 
Q] After the first seven points plotted in Figure Q] the 
F10.7 has been smoothed as per the frequency shifts. 
To determine the scale on which the proxy axis was 
plotted, a linear least-squares fit between the observed 



frequency shifts and the activity proxy was performed. 

The rise and decline of cycle 23 appears to be slow 
in comparison to cycle 22. The oscillation frequencies 
observed in the current unusual solar minimum are sig- 
nificantly lower than the previous minimum in 1996 
and the most recent data indicate that we are still on 
a downward trend. Notice that the shapes of the pre- 
vious and current solar minima are very similar, with 
both appearing to show a double minimum. This im- 
plies that we may still have to wait some time before 
the rising phase of cycle 24 begins. It is also interest- 
ing to note that, if account is taken of the lower quality 
of the data observed before 1990, the minimum be- 
tween cycles 21 and 22 is also deeper than the last so- 
lar minimum. This could be because both the mini- 
mum between cycles 21 and 22 and the current solar 
minimum (between cycles 23 and 24) correspond to 
the same phase of 22-year magnetic Hale cycles. 

It is clearly apparent from Figure [1] that there is an 
unusually large difference between the activity proxy 
and the frequency shifts on the declining phase of cycle 
23. Notably, the frequency shifts are still observed to 
be changing throughout the solar minimum, which is 
in contrast to the activity proxy data that show little 
structure in its variation with time since 2006. 

Figure [2] allows a comparison between the BiSON 
frequency shifts and two proxies of the Sun's activ- 
ity, other than the F10.7. As noted above the ISN is a 
measure of the number of active sunspot groups and 
individual s pots on the visible di sk. The Mgll core-to- 
wing ratio dViereck et al. 2001) is determined from 
space-based UV spectral irradiance measurements. 
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Fig. 2. — Comparison between the frequency shifts observed in the BiSON data and two different activity proxies. 
Top: overplotted in red is the ISN. Bottom: overplotted in red is the Mgll core-to-wing ratio. 



The Mgll H and K line cores originate in the chro- 
mosphere and the ratio of the core size to that of the 
more stable background gives a robust indication of 
the chromospheric activity. 

The F10.7 and the Mgll core-to-wing ratio appear 
to show better agreement with the frequency shifts 
than the ISN. This b ehavior has been observed before 
( Chaplin et alj2007l) and occurs because the frequency 
shifts are sensitive to both the strong and weak com- 
ponents of the Sun's magnetic flux. The F10.7 and the 
Mgll core-to-wing ratio show similar relative sensitiv- 
ity to the different components of the Sun's magnetic 
flux, whereas, the ISN is predominantly sensitive to 
the strong component of the Sun's magnetic field. 

For the purposes of this Letter it is most important 
to notice that there are large discrepancies between, re- 
spectively, all three of the proxies plotted in Figures Q] 
and [2] and the frequency shifts in the declining phase 
of cycle 23. Furthermore, the structure that is observed 
in the frequency shifts in the current solar minimum 
is not replicated by any of the proxies. It is also im- 
portant to note that the seismic minimum is decidedly 
lower than the minimum predicted by the three prox- 



ies. The oscillations are sensitive to the conditions be- 
neath the solar surface, whereas the activity proxy is 
a measure of the magnetic field that is present on the 
surface. We therefore assume that there are processes 
happening in the solar interior, which have yet to man- 
ifest themselves at the surface. 

2.1. Short-Term Variations in the Frequency 
Shifts 

We now turn our attention to the shorter-term vari- 
able structure in the frequency shifts visible on top 
of the general 11 -year trend, which has a period of 
around two years. Similar "quasi-biennial" variation 
of activity proxies has been noted before in data from 
the green cor onal emission line at 530. 3 nm at high 
solar activity dVecchio & Carbonal2008 ). These au- 
thors used variations near the solar equator and poles 
and evidence was found of quasi-biennial variabil- 
ity at both low and high latitudes. An explanation 
of such quasi-biennial behavior has been put forward 
in terms of two dif ferent types of dynamo op erating 
at diffe rent depths ( Benevolonskaya 1998al a). Else- 
where ( Saar & Brandenburg! 2002L it has been pos- 
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sible to identify shorter secondary periods in the ac- 
tivity cycles of some stars and it is conceivable that 
asteroseismology may be able to detect such effects 
in data sets which will become available, for ex- 
ample, from the recently launche d Kepler satellite 
dChristensen-Dalsgaard et al.ll2008l) . 

In spite of the known behavior of the proxies, this 
is the first time that the quasi-biennial variability has 
been noted in seismic data. Furthermore, the signature 
of the two-year signal in the activity index was previ- 
ously restricted to times of moderate to high solar ac- 
tivity. For the first time, we are seeing this short-term 
variability at low solar activity. To investigate this fur- 
ther, a sine wave was fitted to the frequency shifts ob- 
served after 1985 April (see Figure[TJ. The sine wave 
took the form 

S(6t) = A sin(27rA 1 <5f + A 2 ) + A 3 6t + A 4 , (1) 

where 5t is the time in years after 1985 April 22, which 
is the start date of the first 54 d time series. The A, 
coefficients are given in Table \T\ The sine wave was 
calculated using a weighted least-squares fit to the un- 
smoothed frequency shifts. In equationQ] the first col- 
lection of terms accounts for the sine wave structure of 
the frequency shifts. The second and third terms ac- 
count for an offset which linearly decreases with time. 

We have determined the residuals between the best- 
fitting sine wave and the observed frequency shifts 
and the results are plotted in Figure [3] Flux residu- 
als were also determined between a sine wave, which 
was scaled in amplitude using the linear fit between the 
flux and the frequency shifts, and the F10.7. Clearly, 
there is a lot of structure in both sets of observed resid- 
uals and since ~ 2002 the frequency-shift residuals are 
quite substantial in size. Large discrepancies between 
the two sets of residuals are clearly evident. A peri- 
odogram of the data from the last two cycles shows 
significant periods of about two and three years with 
greater than 99% confidence, but the phase of the sig- 



Table 1: The A,- coefficients used in calculating equa- 
tion Q] 
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nal is not locked to that of the underlying 1 1 -year cy- 
cle and several incidences of phase jumps are evident, 
meaning that periods are not easy to determine. The 
three-year periodicity is not visible in a periodogram 
of the activity proxy residuals. 

3. Summary 

The solar-cycle shifts that are observed in p-mode 
frequencies are usually well correlated with proxies of 
the Sun's activity, such as the 10.7 cm radio flux. How- 
ever, in the declining phase of cycle 23 and the current 
solar minimum we find unusually large differences be- 
tween the frequencies observed in BiSON data and the 
activity levels. The current cycle minimum indicated 
by the helioseismic data is significantly deeper than 
the minima observed by the activity proxies and the 
structure that is clearly evident in the frequency shifts 
is not replicated in the proxy data. We also observe a 
quasi-biennial signal in the p-mode frequencies at all 
activity levels. Interestingly, this signal is only visi- 
ble at /i/g/z-activity levels in the proxy data. As the 
frequency shifts respond to conditions beneath the sur- 
face of the Sun whereas the proxies respond to changes 
at or above the surface, we suggest that these differ- 
ences may be caused by changes in the magnetic flux 
that have yet to manifest at the surface. It is also possi- 
ble that the magnetic flux responsible for the discrep- 
ancies between the frequency shifts and the activity 
proxies will never reach the Sun's surface. The anal- 
ysis presented here was based on averages made over 
groups of modes (0 < I < 3, 14 < n < 23). Further 
work on individual modes may allow one to isolate the 
location of the variability because each mode shows 
a different sensitivity to the latitudinal distribution of 
the surface activity. Such an analysis is currently in 
progress. 

In the seismic data, the previous solar minimum 
exhibited a double minimum and it appears that the 
current solar minimum is showing a similar structure. 
Furthermore, the most recent p-mode frequencies indi- 
cate that the current minimum is still declining. There- 
fore, it may still be some time before we observe the 
rising phase of cycle 24. There have been suggestions 
that this recent strange behavior of the Sun is indica- 
tive that the curren t grand maximum is about to end 



( Abr eu et al.l 12008). That would indeed be an occur- 



rence of great significance. Although our results can- 
not predict whether this is true they do indicate that the 
next solar cycle should be observed with great interest. 
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Fig. 3. — Residuals observed between the best-fitting sine wave described by equation [TJ and the observed frequency 
shifts. Overplotted in red, and using the right-hand axis, are the residuals between the sine wave, whose amplitude 
was scaled so as to be appropriate to the flux measurements, and the F10.7. The F10.7 residuals have been scaled using 
a linear fit to the observed frequency shifts. 
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